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ABSTRACT 

We undertake a theoretical study of the near-infrared (NIR) hghtcurves of Type la supernovae 
(SNe la). In these bands, the hghtcurves are distinguished by a secondary maximum occurring roughly 
20 to 30 days after the initial one. Using time-dependent multi-group radiative transfer calculations, 
we calculate the UBVRIJHK-band hghtcurves of model SN la ejecta structures. Our synthetic NIR 
hghtcurves show distinct secondary maxima, and provide favorable fits to observed SNe la. We offer a 
detailed explanation of the origin of the NIR secondary maximum, which is shown to relate directly to 
the ionization evolution of iron group elements in the ejecta. This understanding provides immediate 
intuition into the dependence of the NIR hghtcurves on the physical properties of the ejecta, and in 
particular explains why brighter supernovae have a later and more prominent secondary maximum. 
We demonstrate the dependence of the NIR hghtcurves on the mass of ^^Ni, the degree of ^^Ni mixing, 
the mass of electron caputre elements, the progenitor metallicity, and the abundance of intermediate 
mass elements (especially calcium). The secondary maximum is shown to be a valuable diagnostic 
of these important physical parameters. The models further confirm that SNe la should be excellent 
standard candles in the NIR, with a dispersion ^0.2 mag even when the physical properties of the 
ejecta are varied widely. This study emphasizes the consummate value of NIR observations in probing 
the structure of SNe la and in furthering their cosmological utility. 

Subject headings: radiative transfer - supernovae: general 



1. INTRODUCTION 

The extensive monitoring of Type la supernovae 
(SNe la) has typically focused on optical band obser- 
vations, with only occasional ventures into t he infrared 
l|Kirshner et al.l ll973': "Ehas et alllToSll [l98^. The sit- 
uation is now beginning to change, and near-infrared 
(NIR) observations over the two or three months fol- 
lowing explosion are becoming more common, at least 
for the more nearby events. About a dozen SNe la 
have published, well observed JHK-band hghtcurves 
(e.g., EvIciklc 2000; Hernandez c t al. 20 00; Valentini et aL 
20mt jCandia et alJ 1200^ iKrisciunas et alJ I2001L 12001 
2004blcf) ! NIR spectra have been obtained for sev- 



eral objects as w ell (Mei kle et al.l 199^ Bowers et al, 
[1997':'Rud v et aUlToOa iHofiich et al II2002I: iMa rion et al 
[2003). Future SN la surveys promise to gather NIR 
hghtcurves for a statistically interesting sample of events. 
With the broadening of our wavelength horizons comes 
the potential to probe the nature of SN la explosions in 
new ways. 

Observational studies of SNe la in th e NIR have often 
empha sized their cosmological utility ijKrisciunas et al.l 
l2nn4aD . Normal SNe la are found to be excellent stan- 
dard (as opposed to calibrated) candles in the NIR, with 
an intrinsic dispersion of less than 0.2 mag in the J,H, 
and K bands (Ehas ct al. 1985; Meikle 2000). Moreover, 
dust extinction is smaller by a factor ~ 5 in the NIR as 
compared to the V band, largely eliminating uncertain- 
ties in the reddening corrections. The primary challenge 
for the cosmology studies, naturally, is that SNe la are 
much fainter and harder to observe in the NIR. 

The NIR hghtcurves of SNe la possess a morphology 
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distinct from those at optical wavelengths. The I through 
K-band hghtcurves are distinguished by the presence of 
a secondary maximum, occurring roughly 20 to 30 days 
after the initial one. A corresponding "shoulder" can of- 
ten be seen also in the R and V-band hghtcurves. When 
the light from all bands is suitably integrated, the bolo- 
metric hghtcurves sometimes show an inflection at these 
times as well (Contardo et al. 2000). 

In the I-band, the properties of the secondary max- 
imum are found to correlate with the lightcurve de- 
cline rate (and hence peak luminosity) of the su- 
pernova, being more prominent and occurring later 
in the broa der /brighter SNe la (Hamuv et al. 1996^ 
iNobili et al.l(200lTl) . Very subluminous objects may lack 
a secondary maximum entirely. Similar trends appear 
in the J, H, and K band hghtcur ves, although there are 
interesting exceptions to the rule ijKrisciunas erani20(nt 
^andia ct al. 2003). 

Occasional theoretical stud ies have touched up on the 
NIR hghtcurves of SNe la. iHofiich et all ljl995(l com- 
puted I-band and NIR hghtcurves for several delayed- 
detonation models, many of which displayed secondary 
maxima. The authors explained the double-peaked be- 
havior as a "temperature-radius" effect, in which the ex- 
pansion of the photosphere c ompensates for the de clin- 
ing temperature in the ejecta. iWheeler et alJ l)1998(l and 
IHofiich et al.l l)2002|) have further demonstrated the value 
of NIR spectr a in diagnosing th e phy sical conditions 
in the ejecta. iPinto fc Eastmanl (|2P00) described the 
NIR secondary maximum as the release of pre-existing, 
trapped radiative energy due to a sudden decrease in 
the flux mean opacity. They placed special emphasis 
on the role played by stable iron group elements at 
the eje cta center in this process. Most significantly, 
IPinto fc Eastman (20Q0) recognized the importance of 
the ionization state of the ejecta, in particular the in- 
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creased NIR emissivity of singly as opposed to double 
ionized iron group elements. This last idea is confirmed 
and extended in the models studied here. 

Despite the theoretical insights, the origin of the NIR 
secondary maximum and its dependence on the SN ejecta 
properties have remained in many ways, and to many 
people, obscure. Moreover, theoretical models of SNe la 
have often had difficulty in fitting the distinctive double 
peake d I-band lightcurv e s (e.g ., iHocflich & Khokhlo 
19961 iPinto fc EastmanI I2000I: iBhnnikov fc Sorokin; 
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2004|) whereas very f ew NIR lightcurve c alculations have 



been attempted fc.f. iHoflich et al]l 19951) . This paper fo- 
cuses exclusively on theoretical models of the far-red and 
NIR lightcurves of SNe la. We calculate the synthetic 
lightcurves of parameterized ejecta configuration using 
the time-dependent multi-dimensional radiative transfer 
code SEDDNA Our model lightcurves show distinct 

secondary maxima and provide reasonable fits to the ob- 
served R, I, J and (to a lesser extent) H and K band 
hghtcurves of SNe la (gl). 

Close examination of our model calculations provides 
a clear-cut explanation for the NIR secondary maximum 
in SNe la The double-peaked behavior is related 

directly to the ionization evolution of iron group ele- 
ments in the SN ejecta. In particular, the NIR emis- 
sivity of iron/cobalt increases sharply at a temperature 
T w 7000 K marking the transition between the singly 
and doubly ionized states. The eventual cooling of the 
iron-rich layers of ejecta to this temperature is thus ac- 
companied by a sudden increase in emission at far red 
and NIR wavelengths. 

Understanding the origin of the secondary maximum 
provides immediate intuition into the dependence of the 
NIR lightcurves on the ejecta properties, and allows for 
theoretical explanations for the empirically established 
correlations (S©. In particular we study the effect of 
variations in the mass of ^^Ni, the degree of mixing of 
^^Ni, the mass of electron capture elements, the pro- 
genitor metallicity, and the abundance of intermediate 
mass elements (especially calcium). All effects are con- 
firmed with illustrative theoretical models. The NIR sec- 
ondary maximum is shown to be a valuable diagnostic of 
these important physical parameters. The models fur- 
ther demonstrate that SNe la are indeed excellent stan- 
dard candles in the NIR, even when the physical proper- 
ties of the ejecta are varied widely (IQ. 

2. THE NEAR-IR EMISSION OF SNE lA 

Before turning to the lightcurve models, we consider 
the source of NIR emission in SN la atmospheres. The 
luminosity of a SN la is powered entirely by the decay of 
radioactive isotopes synthesized in the explosion, in par- 
ticular 56Ni in the decay chain ^^Co ^^Fe. 
The energetic products of the decay (primarily gamma- 
rays) heat the SN ejecta, which radiate thermally. The 
bulk of the radiation occurs initially at ultraviolet (UV) 
and blue wavelengths. Direct emission at NIR wave- 
lengths (from bound-free and free-free sources) is typi- 
cally of minor significance. 

A portion of the blue/UV radiation pervading the SN 
ejecta is degraded to red and NIR wavelengths by in- 
teraction with lines. The primary re distribution mech- 
anism is fiuorescence (Hoflich .19951 IPinto fc EastmanI 
120001) whereby a UV photon excites a high-energy atomic 
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Fig. 1. — Opacity as a function of wavelength of a parcel of 
cobalt/iron ricli ejecta gas at density p = 10"^'' cm~^, tempera- 
ture T = 15000 K and time since explosion toxp= 20 days. The 
composition is 75% ^^Co, 24% iron and 1% 56Ni. LTE level popu- 
lations are assumed. The thick black line shows the bound-bound 
expansion opacity (Equation calculated using the Kurucz GDI 
atomic linelist (~ 42 million lines), while the thin red lino is that 
calculated using the Kurucz GDI linelist (~ 500,000 lines). Note 
the large differences in the near-IR. The horizontal dashed line 
marks the electron scattering opacity level. 



transition, then de-excites via a cascade of longer wave- 
length transitions (the reverse process is also possible and 
occurs as well; Mazzali & Lucv 1993). The wavelength 
redistribution of radiation in atomic lines is responsi- 
ble for virtually all of the NIR emission observed from 
SNe la. The high velocity gradient in the ejecta Doppler 
broadens the radiation from each individual transition, 
and the radiation from literally millions of lines blurs to- 
gether to form a pseudo-continuum emission in the NIR. 

We can calculate, in an approximate way, the NIR 
opacity and emissivity due to the large complex of 
lines. A convenient way to capture the average effect 
is through an expansion opac ity formalism (Karo ct aD 
119771 lEastman fc Pintolll99^ . in which individual lines 
are summed together to form a pseudo-continuum opac- 
ity 
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where texp is the time since explosion. The sum runs over 
all lines in the wavelength bin of size AA at wave length A. 
The quantity t; is the line Sobolev optical depth ijMihalad 
B978). ^ 

In Figure ^ we show the expansion opacity calcu- 
lated for a parcel of iron/cobalt rich ejecta gas at den- 
sity p = 10"^'^ cm^-^, temperature T — 15000 K and 
time icxp= 20 days, typical of the inner layers of SNe la. 
We use the Kurucz GDI atomic linelist containing ~ 42 
million lines, and compute the atomic level populations 
assuming local thermodynamic equilibrium (LTE). The 
iron group expansion opacity increases sharply towards 
the blue, due to the much larger number of iron group 
lines occurring at shorter wavelengths. 

In Figure ^ we also compare the expansion opac- 
ity computed using the much smaller (but commonly 
applied) Kurucz CD 23 linelist, containing only ~ 
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500, 000 lines. Although the two lists give similar re- 
sults in the blue, order of magnitude differences appear 
in the NIR. This emphasizes the large number (~ 5 mil- 
lion) of relatively weak lines that may contribute to the 
NIR emission of SNc la. 

Assuming LTE, the monochromatic emissivity r]\ is 
given by the opacity acxp times the blackbody function 
B\{T). Because we are ultimately interested in the to- 
tal NIR emission over a broadband wavelength region, we 
define an average emissivity per gram (units ergs/s/g/A) 
over the broadband filter F by convolving the monochro- 
matic emissivity with the appropriate filter transmission 
function 0_f(A) 



Vf{T) = — ^ 



P 



(2) 



Although this expression formally represents a thermal 
emission source, under certain conditions it well approx- 
imates the operative line fluorescence process as well. In 
the limit that the time-scales for radiative transitions 
are short compared other time-scales of interest, the ra- 
diation field in the lines approaches the thermodynamic 
limit embodied by the Planck function. These are the 
conditions that drive the line source function and level 
populations close to that of LTE, and should hold ap- 
proximately at least for the c omplex iron group species 
in the deeper layers of ejecta ijBaron et alJll996|) . Given 
the dominance of the radiative rates, the temperature T 
appearing in Equation [21 should be considered more rep- 
resentative of the radiation field, rather than the electron 
gas. However, the weak coupling through bound-free and 
free-free processes may bring the gas and ra diation into 
equilibrium as well (Pinto & Eastmanl l2000(l . 

Note that, even assuming LTE, the wavelength depen- 
dence of the emissivity is not a smooth blackbody func- 
tion, but depends upon the the number of lines with 
substantial opacity in a given part of the spectrum. As a 
result, the NIR emissivity shows a very interesting non- 
linear dependence upon temperature. In Figure 12 we 
see that the mean I-band emissivity (77/) of iron/cobalt- 
rich gas peaks sharply at temperatures where the gas is 
near an ionization edge. The most prominent peak is 
at a temperature T21 7000 K, marking the transition 
between the singly and doubly ionized states (hereafter 
termed the 2 — > 1 ionization edge or front). 

The interesting dependence of the NIR emissivity on 
temperature has a simple explanation in terms of the 
atomic physics. NIR radiation arises from transitions be- 
tween relatively closely spaced atomic energy levels. The 
atomic levels are generally more closely spaced at high 
excitation energy (i.e. high atomic number n). Thus, 
for a given ionization state, the NIR emission increases 
with temperature as the excited levels become increas- 
ingly populated. When the temperature becomes hot 
enough to ionize the species, however, the NIR emissivity 
drops rapidly. The ionized species has a higher overall 
energy scale (being more tightly bound) and hence its 
excited levels are not highly populated. If the tempera- 
ture is increased further, this behavior is repeated until 
the next ionization stage is reached. We emphasize that 
this emissivity dependence has nothing to do with the 
energy released from recombination itself, which makes 
a negligible contribution to the overall energy budget. 
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Fig. 2. — Mean I-band emissivity (as a function of temper- 
ature) for the iron/cobalt gas discussed in Figure at density 
p = 10"^'' cm"'' and time since explosion tcxp= 40 days. Cal- 
culations use Equation |^ the Kurucz GDI atomic linelist, and 
the Bessell 1 1990) I-band filter profile. The shading illustrates the 
ionization fraction of the gas, which changes very sharply with 
temperature at the edges marked in the figure. A strong peak in 
emissivity is seen at the transition between the singly and doubly 
ionized states (2 1) occurring at T ss 7000 K. 



For completeness. Figure |31 shows the emissivity for 
several wavelength bands and for compositions of pure 
iron, pure cobalt, and a mixture of silicon and sulfur. 
The emissivity of the iron-group species is qualitatively 
similar in all NIR bands. A peak at T21 = 7000 K oc- 
curs also in the optical B and V bands, although it is 
much less prominent than in the NIR. The relative size of 
the emissivity peaks depends upon the details of atomic 
structure, and thus varies somewhat between iron and 
cobalt and in the different wavelength regions. In gen- 
eral, the iron-group emissivity exceeds that of interme- 
diate mass elements such as silicon and sulfur, except 
at very high temperatures. This is because the complex 
iron group valence electron structure leads to a extremely 
large number of line transitions. 

The emissivity behavior shown in Figure |2] turns out 
to be the key to understanding the NIR secondary max- 
imum in SNe la. When the iron rich layers of SN ejecta 
begin cooling to T21 « 7000 (which occurs roughly 35 
to 40 days after the explosion) the gas becomes very ef- 
fective in redistributing the blue/UV radiation to longer 
wavelengths. This leads to the rebrightening of the NIR 
lightcurves. During these phases, one expects the NIR 
emission to be dominated by the relatively thin shell of 
material located near the ionization front at T21. As the 
ejecta progressively cools, this shell recedes deeper into 
the ejecta, following the inward propagation of the 2^1 
recombination. 

The secondary maximum in the NIR lightcurves is thus 
identified with the onset and recession of the 2 — > 1 ion- 
ization front into iron-rich ejecta. In the next few sec- 
tions we delineate the process in detail, and show how it 
serves as an interesting probe of the SN ejecta structure. 

3. THEORETICAL SNE lA LIGHTCURVE MODELS 

A short time after the explosion of a SN la (~ 1 
minute) hydrodynamic and nucleosynthetic processes 
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Fig. 3. — A more complete plot of the mean NIR emissivity 
demonstrated in Figure|21 this time shown for each of the BVIJHK 
bands and for three compositions: pure iron (solid black), pure 
cobalt (dashed red) and a mix of 70% silicon and 30% sulfur (dot- 
dashed blue). The density is p = 10"^"* cm~^ and the time since 
explosion tcxp= 40 days. 



abate and the ejected material reaches a phase of nearly 
free expansion. Thereafter, the essential theoretical chal- 
lenge becomes to simulate the production and diffusion 
of photons through the hot and optically thick ejecta 
- i.e., the radiative transfer problem. Here we calcu- 
late broadband model lightcurves by solving the trans- 
fer equation using the multi-dimensional time-dependent 
Monte Carlo code SEDDNA (Kasen et al. 20Q<Q. 

3.1. Model Ejecta Structures 

To demonstrate the basic dependencies of the NIR 
lightcurves, we calculate synthetic SEDDNA lightcurves for 
simple, hand-constructed SN la ejecta models. In the 
free-expansion phase, the velocity of the ejecta is homol- 
ogous and everywhere proportional to radius: r — vt^xp 
where icxp is the time since explosion. Given the self- 
similar nature of the flow, velocity is used as the spatial 
coordinate in the simulation. 

Polarization observations of SNe la demonstrate that 
the as ymmetry of the ejecta is typically not large 
ijWang et al. 2003). We therefore only consider spheri- 
cally symmetric ejecta structures in this paper. Devia- 
tions from sphericity will not effect the basic radiative 
transfer effects identified here, though they will lead to 



some orientation dependence of the model lightcurves. 

In 1-dimension (1-D) hvdrodvnamic explos ion models 
of SNe la such as W7 l|Nomoto et al.lll984D . the final 
ejecta density structure is well characterized by an ex- 
ponential. We therefore take the density profile of our 
ejecta models to be 



M, 



ch 



e-Kp{-v/Ve), 



(3) 



where Mch is the Chandrasekhar mass and the velocity 
scale Ve — 2750 km s~^ is determined by setting the ki- 
netic energy KE = GMchV^ equal to the energy released 
in burning 1.25 Mq of the original carbon/oxygen white 
dwarf (viz. 1.6 x 10^^ ergs). 

In analogy to well-known 1-D explosion models such as 
W7, we tailor a four-zone stratified compositional struc- 
ture. For our fiducial model these zones are, from the 
center out: (1) 0.05 Mq of stable iron peak elements 
(specifically ^''Fe and ^^Ni, typical of electron capture); 
(2) 0.6 Mq of radioactive ^^Ni; (3) 0.6 Mq of intermedi- 
ate mass elements (IME) specifically. Si, S, Ar and Ca, 
in proportion to their solar abundances (4) 0.15 Mq of 
unburned carbon-oxygen with solar metallicity. 

To smoothly connect these zones, we use a Gaussian 
function that allows us to control the degree of mixing. 
For example, mixing between the ^^Ni and IME layers 
is assumed to occur over a mass range defined by the 
free parameter am ■ Above the interface between the two 
zones (mass coordinate m„ = Mpc + Mni — aNi/2) the 
nickel abundance is taken to fall off as 



^Ni = exp(-(M - m„)^/aNi), 



(4) 



while the IME abundance is Xime = 1 — -''^Ni- Similar 
mixing parameters are defined for the iron and ^^Ni inter- 
face (ape) and the IME and unburned interface (oime)- 
In our fiducial model we use ajsn— 0.1, apc— 0.05, and 
aiME= 0.1. In fact, only am is of real consequence to the 
NIR lightcurves ( ^^ . 

The rough compositional structure described, while 
only suggestive of detailed explosion models, is quite ad- 
equate for our purposes. The NIR lightcurves are not 
sensitive to the detailed trace abundances, with one im- 
portant exception - as we discuss in ii6.4l small amounts 
of calcium can potentially affect the I-band, due to the 
overwhelming strength of the Ca II IR triplet feature. In 
all other bands, the overriding factor is the amount and 
distribution of iron- group elements. 

3.2. Technical Considerations of the Radiative Transfer 

Given a homologously expanding SN ejecta structure, 
the time-dependent SEDONA code calculates high resolu- 
tion synthetic spectra at each day, from day one to sev- 
eral months after the explosion. Broadband lightcurves 
are then constructed by convolving the spectrum at each 
time with the appropriate filter tra nsmission functions. 
We use the filter profil es of iBesselll il990) for the opti- 
cal bands and those of .Persson et al.l (^1998) for the JHK 
bands. 

SEDONA includes a detailed treatment of gamma-ray 
transfer to determine the instantaneous energy deposi- 
tion rate from radioactive ^^Ni and ^®Co decay. Radia- 
tive heating and cooling rates are evaluated from Monte 
Carlo estimators, and the temperature structure of the 
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ejecta determined by iterating the model to thermal equi- 
librium. The present calculations use 120 radial zones, 
100 time points, and 5000 wavelength groups. Resolu- 
tion tests confirm the adequacy of this gridding for the 
problem at hand. 

Several significant approximations are made in SEDDNA, 
notably the employment of LTE level populations and 
ionization. In addition, bound-bound line transitions are 
treated using the expansion opacity formalism (Equa- 
tion^ and an approximate two- level atom approach to 
wavelength redistribution. The ratio of absorption to 
scattering in lines is calibrated from detailed atomic mod- 
els and is always close to one. Special care, however, is 
taken for the calcium lines, which are assumed to be 
pure scattering for the reasons discussed in ^6A\ Note 
that SEDONA allows for a direct Monte Carlo treatment 
of line fluorescence, but due to computation al constraints 
this fu nctionality is not exploited here. See lKasen et al.l 
1)20061^ for further code description and verification. 

In the tenuous atmospheres of SNe la, the microscopic 
conditions assuring the establishment of LTE level pop- 
ulations and ionization (namely, the dominance of col- 
lisonal rates) are in fact not met. A direct treatment of 
the non-LTE physics is thus highly desirable. Unfortu- 
nately, a solution of the non-LTE rate equations includ- 
ing the > 5 million potentially important line transitions 
challenges the computational capacity of even the most 
advanced time-independent spectrum codes. In a fully 
time-dependent light curve calculation, the adoption of 
LTE at some level appears inevitable. Fortunately, a 
wide range of theoretical calculations confirm the ade- 
quacy of LTE in reproducing the essential spectral and 
photometric prop erties of SN la models in both the opti- 
cal and NIR, (e.g. iHoflich et aLllTool iBaron et al.lll99fil: 
iPinto fc Eastmanll2000HWheeler et al.l ll998^. In partic- 
ular, the crucial iron-group species responsible for con- 
tinuum formation have extremely complex atomic struc- 
tures and numerous closely spaced atomic levels, and 
may therefore roughly approximate their equilibrium dis- 
tributions. However, in other atomic species responsible 
for individual strong line features non-LTE effects may 
be significant (especially calcium, see ii6.4|l . In addition, 
because the ionization state of the ejecta is of particu- 
lar importance to the NIR lightcurves, LTE becomes in- 
creasingly suspect at later times (toxp ^ 70 days), when 
the non-thermal ionization by radioactive gamma-rays 
becomes significant (SQ. 

An additional, perhaps more serious issue for the 
transfer calculations is the (in)adequacy of the available 
atomic data. Because NIR emission occurs in millions of 
relatively weak atomic lines, using a complete and accu- 
rate atomic linelist is critical. We have f o und t he com- 
monly applied CD 23 linelist of iKuruc 3 lll99l (which 
contains nearly 500,000 lines) to be insufficient for calcu- 
lation of the IJHK-band lightcurves. Upgrading to the 
Kurucz CD 1 linelist (with ^ 42 million lines), affords 
significant improvement, however even this list is likely 
incomplete and/or inaccurate to some degree. The inad- 
equacy of the line data at long wavelengths likely plagues 
our synthetic NIR lightcurves, as we discuss further in 21 

4. MODEL NIR LIGHTCURVE OF A NORMAL SNE lA 

We have calculated the full spectral evolution and 
broadband lightcurves for our fiducial SN la model con- 



taining 0.6 M© of ^^Ni. FigureQlcompares the model op- 
tical and NIR lightcurves to observations of SN 2001el, 
a photometrically normal S Ne la with extensive cover- 
age obtained bv Krisciunas et al.l l)2003(l . The observed 
lightcurves have been corre cted for dust exti nction us- 
ing the estim ates A„ = 0.57 llKrisciunas et al.ll2Q031 and 
= 2.88 ijWang et all 12003(1 . We select a distance 
modulus to SN 2001el of ^ = 31.4 in order to align 
the peak of the observed y-band lightcurve with that 
of the model. This improves t he visual com parison of 
the lightcurve shapes. In fact, iKrisciunas et al. (2003) 
estimate a slightly smaller distan ce of 31.29 ± 0.08 us - 
ing the width- luminosity relation o flPhillips et al.ll)1999l) . 
This suggests that our synthetic lightcurves are perhaps 
~ 0.1 mag brighter than the observed, and that a model 
with slightly lower A'/ni (e.g., 0.55 Mq) may be more 
appropriate for SN 2001el. 

The characteristic double-peaked morphology is clearly 
apparent in all our model NIR lightcurves, with a dis- 
tinct and prominent secondary maximum occurring in 
the I,J,II, and K-bands. A corresponding "shoulder" is 
also seen in the R band, and in the bolometric lightcurve 
as well (Figure El). On the whole, the model V,R,I and 
J-band lightcurves do a favorable job of fitting those of 
SN 2001el, at least for epochs tcxp < 60 days. The 
models get progressively worse at later times, but this 
is also when the transfer calculations are expected to be 
less reliable. In the model I-band lightcurve, the decline 
from the first maximum is somewhat too rapid and the 
secondary maximum occurs slightly earlier than in the 
observations. We note that the prominence and tim- 
ing of the secondary maximum is quite sensitive to the 
model parameters (see and also show s wide varia- 
tion in the observations l)Nobili et al.ll2005l) . A thorough 
exploration of the model parameter space would likely 
yield an improved fit to SN 2001el, perhaps in the pro- 
cess illuminating its particular physical characteristics. 
In addition, polarization observations of SN 2001el indi- 
cate th at the ejecta. was not perfectly spher ically sym- 
metric fWang et all 120031: iKasen et alJl2003H . therefore 
aspherical ejecta structures may need to be considered 
to exactly account for all of the lightcurve properties. 

The model does a poorer job in the H and K bands. 
In the observations, the H- and K-band lightcurves are 
relatively flat for roughly twenty days after B-band max- 
imum (icxp= 18 days), whereas the models, in contrast, 
show a strong rise to the secondary maximum during 
these phases. This is because the model underestimates 
the luminosity of the initial maximum relative to the sec- 
ondary one. Attempts to improve the H and K fits by 
varying the model ejecta structure proved unsuccessful. 

It is possible that the relatively poor H and K band 
lightcurve fits relate to an inadequacy of the LTE ap- 
proximation, or signal a missing ingredient in our model 
ejecta structures. However, the primary reason is likely 
the inadequacy of the atomic linelist. We demonstrate 
the importance of the atomic data by recomputing the 
lightcurves using the restricted linelist Kurucz CD23 
(containing only 500,000 fines) and comparing to those 
of our standard linelist (Kurucz CDl, with 42 million 
lines). The differences in the IJHK lightcurves are con- 
siderable (Figure El, especially at early times. This is 
not surprising, given the large differences noticed in the 
opacity calculations of Figure H 
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On the whole, the CD23 IJHK hghtcurves do a dra- 
matically poorer job of fitting the observations. In par- 
ticular, the initial maximum is much to weak in all bands. 
The problem can be traced directly to the relative lack 
of higher ionization, longer wavelength iron group lines 
in the CD23 list. The use of an inadequate linelist 
may be one reason why previous theoretical models have 
sometimes had difficulty fitting the I-band lightcurves of 
SNe la. Although GDI offers significant improvement 
over CD23, it is natural to suspect that were an even 
more extensive linelist available, our model H and K 
band lightcurves might show significant improvement. 

Although not the focus of this paper, we also present 
NIR spectra of the fiducial model at a few different 
epochs in Figure [3 Near maximum light, the NIR spec- 
trum is quite featureless, with the occurrence of only a 
few weak line features from intermediate mass elements. 
At later times, prominent emission features from blended 
iron-group lines arise near 1.6 and 1.8 /j,. The spectra in 
Figure 13 show the same general features identified in the 



model calculations of I Wheeler et al.l l)1998t) . The reader 
can look there for a detailed spectroscopic analysis and 
comparison with observations. 

5. EXPLANATION OF THE SECONDARY MAXIMUM 

In our models, we trace the origin of the secondary 
maximum directly to the ionization evolution of iron 
group elements in the ejecta. In particular, as discussed 
in the NIR emissivity increases sharply at tempera- 
ture T21 ~ 7000 K in iron/cobalt gas, marking the tran- 
sition between the doubly and singly ionized states (i.e., 
the 2^1 ionization edge; see Figure For reasons al- 
ready discussed, the iron-rich gas becomes exceptionally 
phosphorescent at this temperature, and is very effec- 
tive in redistributing the pervading UV/blue radiation 
to longer wavelengths. 

Because the temperature in the ejecta decreases with 
radius, the region near the 2 — > 1 ionization front com- 
prises a relative thin shell of material. During the hot, 
early epochs, this region lies in the far outer layers of 
ejecta, but as the SN cools overall the ionization front 
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Fig. 5. — Bolometric lightcurve (solid line) and radioactive en- 
ergy deposition rate (dashed line) for the fiducial model. An in- 
flection in the bolometric lightcurve (corresponding to the NIR 
secondary maximum) is seen beginning at day 35 after explosion. 
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Fig. 6. — Comparison of the IJHK lightcurves of the fiducial 
model calculated using two different atomic linelists: the 500,000 
line Kurucz CD 23 (dashed red lines) and the 42 million line Kurucz 
CD f (solid black lines). Due to the lesser number of lines, the first 
maximum of the CD23 lightcurves is depressed relative to CDl in 
all bands. 



recedes deeper (in a Lagrangian sense) into the ejecta. 
The resulting inward propagating "recombination wave" 
resembles that often described in models of hydrogen- 
rich Type IIP supernovae. However, in contrast to the 
Type HP's, the flux mean opacity in SNe la does not 
change substantially over the ionization front and the 
temperature gradient remains gradual. It is the onset of 
this recombination wave into the iron-rich core of ejecta 
that leads to a sudden increase in NIR emission and the 
re-brightening of the SN in a secondary maximum. 
To delineate this process in detail, we examine the evo- 




1.0x10' 



1.5x10' 2.0x10' 
Wovelength (Angstroems) 



2.5x10' 



Fig. 7. — Near-infrared spectra of the fiducial SN la model at 
three different epochs. Labels indicate time since explosion. 



lution of the NIR emission over time. Figure|Hlis a "comic 
strip" of the spatial distribution of NIR emissivity (in this 
case, J-band) for the fiducial SN la model discussed in 
the last section. Each panel shows the radial contribu- 
tion to the J-band emission at a given time, specifically 
dip{r) = rij{r)r^ dr where r is the radius and rjj is the 
broadband emissivity calculated using Equation |21 This 
figure is meant mainly as a heuristic - the observed lu- 
minosity in the optically thick phases depends not only 
on the instantaneous emission but also the transfer of 
radiation through the ejecta. 

Given the self-similar nature of the flow, velocity is 
used as the radial coordinate in Figure |H1 The red 
hatched region identifies the iron-rich core, which ex- 
tends to roughly 7000 km s~^ in the model. The green 
shaded region marks the region near the 2 — > 1 ioniza- 
tion front (specifically, ionization fraction between 1.5 
and 1.9). The width of this region, which depends upon 
the temperature gradient, is typically a few thousand 
km s""'^. With time, the ejecta cool and the green shaded 
region can be seen moving inward, marking the propaga- 
tion of the 2 1 recombination wave. 

We offer a step-by-step discussion of Figure |H1 Dur- 
ing the early epochs (day 15 after explosion), most of 
the J-band emission comes from throughout the hot iron 
core, which has ionization fraction between 3 and 4. Note 
that the emission is always greatest near the ionization 
fronts, which are noticed as the emission peaks in the fig- 
ure. At these early times, the ejecta is optically thick at 
NIR wavelengths to both electron scattering and bound- 
bound expansion opacity. The rise time to the first J- 
band maximum at texp= 19 days is th us related to the 
diffusion time of photons ljArnettlll982|) . 

Following the first J-band maximum, the hot iron core 
is in the process of recombining from triply to doubly ion- 
ized. At these times (days 20 and 25), the NIR emission 
is dominated by the rather thin shell of material near the 
3 — > 2 ionization front. Meanwhile, the 2 — > 1 ionization 
front forms well outside the iron core and contributes 
little NIR emission. 

By day 30 after explosion, the iron core has nearly 
completely recombined to doubly ionized and, with the 
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Fig. 8. — Comic strip explaining the origin of the secondary maximum in the NIR lightcurves of SNe la. The figure shows the time 
evolution of the radial distribution of the instantaneous mean J-band emission (solid lines). For each panel, the x-axis is the ejecta velocity 
coordinate ranging from zero to 20,000 km 8"^ The red hatched region marks the iron core, which extends to v = 7000 km s ^. The 
green shaded region marks the region near the 2 — > 1 ionization front. A large spike in J-band emission occurs when that front reaches the 
iron rich layers of ejecta, leading to the NIR secondary maximum. 



passing of the 3^2 ionization front, the NIR emissivity 
has decreased significantly. This is the local minimum 
in the J-band lightcurve. The lightcurve would continue 
to decrease monotonically if it weren't for the impending 
onset of the 2^1 recombination wave. Already this 
wave has reached the layers just above the iron core, 
where a small emission peak is visible. 

At day 35 the recombination wave reaches the edge of 
the iron core, and the NIR emission begins to increase 
dramatically. This is exactly the time at which the model 
J-band lightcurve begins to re-brighten. The NIR emis- 
sion is now dominated by the "fluorescent shell" at the 
2 — > 1 ionization front. The ejecta are largely transparent 
in the NIR and radiation emitted from the shell escapes 
almost immediately. Over the next two weeks, the fluo- 
rescent shell recedes further into the iron core, until, at 
day 50, it is completely contained in the iron rich region. 
This marks the peak of the J-band secondary maximum. 

At later times, as the fluorescent shell recedes yet 
deeper in the ejecta, its volume decreases and so does 
the J-band brightness. By day 60, the 2 — > 1 recombi- 



nation wave has completely passed, and the iron core is 
nearly entirely singly ionized. The NIR emission is low 
and continues to decline. 

Why is a secondary maximum seen strongly in the 
NIR, but not in the optical lightcurves? During the 
epochs of the secondary maximum, the ejecta remain 
optically thick at bluer wavelengths (A < 5000 A). The 
onset of the 2 — > 1 recombination of iron group elements 
further increases the opacity in the optical due to the 
development of broad blends of Fe II and Co II lines in 
B and V-bands. Thus optical radiation at the 2 — > 1 ion- 
ization front can not immediately escape the ejecta, but 
continues to diffuse out gradually or to fluoresce to longer 
wavelengths. Meanwhile, the ejecta are transparent in 
the NIR, and photons fluorescing to these wavelengths 
escape straightaway. The net effect of the recombination 
wave is therefore to redistribute flux from the optical to 
longer wavelengths. 

In Figure [S] we see that at the onset of the sec- 
ondary maximum (day 35 after explosion) a slight 
bump also occurs in the model bolometric lightcurve. 
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This inflecti on has been n oted 
hefhtcurve s llContardo et al 



in observed SNe la 
200f][) and has been ex- 



plained bv lPinto &: Eastmar ( 200C ) as the release of pre- 
existing, trapped radiative energy. The ejecta are op- 
tically thick (at bluer wavelengths) for the first couple 
months after explosion, and the average diffusion time of 
photons is a significant fraction of the expansion time. 
A store of trapped radiation continually accumulates in 
the supernova ejecta. It is the gradual leaking of this en- 
ergy reservoir that causes the bolometric luminosity to 
exceed the instantaneous radioactive energy deposition 
rate from day 18 until full transparency is reached near 
day 70. At day 35, the onset of the 2 — > 1 ionization front 
in the iron core suddenly enhances photon escape via flu- 
orescence into the NIR. This opens up a new means of 
release of the trapped SN radiation, causing a slight rise 
in the bolometric lightcurve. 

As an interesting aside, notice that the NIR emissivity 
also peaks when the gas recombines from singly ionized 
to neutral (Figure In principle, the propagation of 
another recombination wave (1 —> 0), could lead to a 
third bump in the NIR lightcurves. In the current mod- 
els, a third maximum is in fact seen at icxp ~ 100 days 
(~ 80 days after B-band maximum), most prominently 
in the J band. However, the inadequacies of the LTE ap- 
proximation at late epochs make this behavior suspect. 
In particular, the non-thermal ionization from radioac- 
tive gamma-ra ys becomes significant at low temperatures 
ijSwartd fl991() and will keep the ejecta from ever go- 
ing completely neutral. This may minimize or delay the 
putative third maximum. In addition, emission by for- 
bidden transitions (not included here) may also become 
significant at these epochs. Meanwhile, most published 
SNe la observations lack NIR data at times > 60 days 
after maximum. 

6. DEPENDENCE ON PHYSICAL PARAMETERS 

The above explanation of the NIR secondary maxi- 
mum provides immediate intuition into the dependence 
of the NIR lightcurves on the physical properties of the 
ejecta. This allows for simple explanations of the ob- 
served trends concerning the timing and prominence of 
the secondary maximum. 

First, the rise to a secondary maximum occurs when 
the 2 — > 1 ionization front first reaches the ejecta layers 
rich in iron group elements. Thus the timing of the sec- 
ondary maximum will depend upon both the size of the 
iron core (or more precisely, its velocity extent), and the 
overall temperature scale of the ejecta. All other things 
being equal, SNe with larger iron cores can be expected 
to have earlier secondary maxima. Hotter SNe can be 
expected to have later ones. 

Second, the luminosity of the second maximum de- 
pends fundamentally upon the size and luminosity of the 
fluorescent shell formed when the 2^1 recombination 
wave has receded into the bulk of the iron core. All other 
things being equal, SNe la with larger iron cores can be 
expected to have a brighter secondary maximum. 

Below we discuss these behaviors in the context of 
SNe la ejecta properties. In particular, we consider the 
effects of the degree of mixing, the ^^Ni mass, the mass 
of electron capture elements, the progenitor metallicity, 
and the mass of intermediate elements produced, demon- 
strating each of these with the models. 
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Fig. 9. — Effect of the outward mixing of ^^Ni on the NIR 
lightcurves. The figure compares the well-stratified fiducial model 
(mixing parameter 0^;= 0.2, solid black lines) to a model with 
enhanced ^®Ni mixing (aNi= 0-6; red dashed lines) and a model 
with a completely homogenized composition structure (blue dot- 
dashed lines). 



6.1. Mixing 

The flducial SN la model studied in the last section 
was compositionally stratified, with iron group elements 
(primarily ^®Ni) concentrated in the central core. The 
mixing of ^^Ni outward into the region of IME effectively 
increases the velocity extent of the iron core. Naturally, 
we expect this to hasten the occurrence of the NIR sec- 
ondary maximum. 

We demonstrate this effect by comparing the NIR 
lightcurves of our fiducial model (mixing parameter 
ONi = 0.2) to one in which ^^Ni is more thoroughly mixed 
outward (oni — 0.5). The result of the mixing (Figure^ 
is to advance the NIR secondary maximum by about 5 
days and to decrease its contrast with the first maximum. 
In the I and J-band lightcurves, the first and second max- 
ima are blended, but still distinguishable. In the H and 
K band lightcurves the double-peaked structure is nearly 
lost entirely. 

In Figure|5|we show also the lightcurves of a fully mixed 
model (i.e., the fiducial model with a completely homog- 
enized compositional structure). In this case, the first 
and second maxima are indistinguishable in every band 
(except for perhaps a small inflection at day 35 in J). We 
conclude that the double-peaked structure observed in 
NIR lightcurves is a direct indication of the abundance 
stratification in SNe la, in particular the concentration 
of iron-peak elements in the central regions. Thus NIR 
observations should be useful in constraining the exact 
degree of mixing in SNe la. 

6.2. Mass o/^^Ni 

The mass of ^^Ni (Mni) powers the luminosity of 
SNe la and is the primary determinate of the peak bright- 
ness of the event. Mni also affects, in part, the B-band 
lightcurve decline rate (or stretch) and is commonly iden- 
tified as the fundamental parameter controlling the well- 
known width-luminosity relation. 

In observations, the properties of the I-band secondary 
maximum are found to correlate with stretch (and hence 
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Fig. 10.— Effect of the mass of ^^Ni (Mni) on the NIR 
hghtcurves. The model Ughtcurves demonstrate variations in Mni 
from 0.1 to 0.9 Mq in 0.1 Mq increments. The less luminous 
models are those with lower Mni- 
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Fig. 11. — Prominence and timi ng o f the secondary maximum 
versus Mni for the models of Figure IkH in the I-band (left) and J- 
band (right). Top panels: Time (in days) of secondary maximum 
minus time of first maximum. Middle panels: Number of mag- 
nitudes that the secondary maximum is brighter than the local 
minimum that precedes it. Bottom panels: Number of magnitudes 
the first maximum is brighter than the secondary maximum. 



peak luminosity), being more prominent and occurring 
later in the broader/brighter SNe la. The origin of these 
empirical relations can be addressed within the theoret- 
ical paradigm we have described. To demonstrate, we 
construct a series of models (based upon the fiducial 
model) in which we vary Mni between 0.1 to 0.9 Mq in 
increments of 0.1 Mq. The sum of Mni and the mass of 
IME is held fixed at 1.25 Mq, such that all models have 
roughly equal kinetic energy. All other model parameters 
have been held fixed. 

As in observations, the model secondary maxima are 
more pronounced in the higher models, disappearing 
altogether for the lowest M^i models (Figure^J. This is 
due primarily to the larger size of the iron core in higher 
Mni models, leading to a larger and more luminous fluo- 



rescent shell at the 2^1 ionization front. The observed 
behavior may therefore be taken as further evidence that 
more slowly declining SNe la have a larger production of 
iron group elements. 

Mni produces two opposing effects on the timing of 
the secondary maximum. A larger Mni leads to higher 
overall temperatures and thus a later development of the 
2 — > 1 recombination wave. This tends to delay the 
secondary maximum. However, a larger Mn! also im- 
plies a larger iron core, which would tend hasten the sec- 
ondary maximum. From Figure [TUl we see that the first 
of these effects dominates, at least for Mni in the range 
0.1 — 0.8 Mq. Thus our models recover the observed tim- 
ing trend, which should be understood fundamentally as 
the physical correlation between the supernova's lumi- 
nosity and its ionization state. 

For the models with low Mni, the secondary maximum 
occurs at early times and begins merging with the initial 
maximum. In our most subluminous example (AfiMi= 
0.1 Mq) the ejecta are so cool that the 2^1 recom- 
bination wave reaches the iron core at texp ~ 20 days. 
For these objects, the secondary maximum is absent, or 
rather it is coincident with the first. 

For reference. Figure 1111 quantifies the timing and 
prominence relations in the I and J bands. We note that 
the observed trends depend to some extent on how the 
quantities are empirically defined. In the J-band, for ex- 
ample, the peak magnitude of the secondary maximum 
increases with Mni relative to the magnitude at the local 
minimum, but decreases relative to the first maximum. 
This refiects complexities glossed over in our discussion 
(such as the parameters affecting the initial maximum) 
and demonstrates the need to compare observations and 
models in a uniform way. Also keep in mind that addi- 
tional parameters (described below) can also have pro- 
found effects on the secondary maximum, and may play 
a role in these relations as well. 

6.3. Electron Capture Elements and Progenitor 
Metallicity 

In addition to ^^Ni, SNe la may synthesize substantial 
amounts of stable iron group elements, in particular ^^Fe 
and ^^Ni. Although this material does not contribute to 
the overall luminosity of the SN, it does modify the radial 
extent of the iron core. Thus we anticipate it having 
important effects on the timing and prominence of the 
secondary maximum. 

Stable iron group material can be produced by two dis- 
tinct processes in SN la explosions: (1) In standard 1-D 
explosion models, nuclear burning at the highest densi- 
ties is subject to electron capture, leading to the syn- 
thesis of 0.5-2.0 Mq of neutronized sta ble iron group 
species at the very center of the ejecta l|Nomoto et alJ 
119841 iThielemann et al]ll986|) . There is some observa- 
tional evidence suggesting such a pure stable iron center 
does exist in SNe la (Hoflich et al. 2004). (2) While burn- 
ing to nuclear statistical equilibrium produces primarily 
^^Ni, some small fraction of stable iron group may also be 
produced, depen ding upon the metallicity of the progen- 
itor w hite dwarf l)Thielemann et aTll 19861 iTimmes et alJ 
I2003|) . For solar metallicity, this fraction would be ~ 5%, 
while for metallicity three times solar the fraction may 
be as high as ~ 25%. The stable iron group material 
arising from the progenitor metallicity would be evenly 
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Fig. 12. — Effect of electron capture elements on the NIR 
lightcurves. The model lightcurves demonstrate variations in the 
mass of stable iron group elements (Mpc) produced by electron 
capture at the ejecta center. The models are Mp(,= 0.0 Mq (black 
line), Mpc= 0.1 Mq (red line), Mpc= 0.2 Mq (blue Une) and Mpo= 
0.3 Mq (green line). 




Fig. 13. — Effect of progenitor metallicity on the NIR lightcurves. 
The black line is the fiducial model with zero stable iron pro- 
duction in the ^^Ni zone (corresponding to near zero metallicity). 
The dashed red line is the model with 25% stable iron production 
throughout the ^^Ni zone (corresponding to a metallicity near 3 
times solar). 



distributed in low abundance throughout the ^^Ni zone. 

We demonstrate each of these effects in turn. First, 
we vary the mass of stable iron group species {Mp^) 
produced by electron capture at the ejecta center. We 
explore values of Mpe from to 3.0 M© while keep- 
ing fixed Mni — 0.6 Mq and the total burned mass 
Mfo + Mm + Mime = 1.25 Mq. Figure [HI shows the 
resulting NIR lightcurves. Because the mass of ^^Ni is 
held fixed, the temperature and ionization evolution is 
very similar in all these models. However, in models 
with larger values of Mpc (and hence larger iron cores) 
the 2 — > 1 recombination wave encounters the layers of 
iron rich material at relatively earlier times. The sec- 
ondary maximum thus occurs as many as ten days earlier 
in models with larger Afpc- In addition, the peak mag- 
nitude of the secondary maximum increases with Afpe, 



a gain due to the larger iro n core size. 

iPinto fc EastmanI l|200nD have suggested that the sta- 
ble iron group material at the ejecta center plays a funda- 
mental role in producing the NIR secondary maximum. 
We note that this is not strictly the case, as a clear sec- 
ondary maximum occurs even when no stable iron is in- 
cluded in the model. 

One finds a similar dependence on the stable iron pro- 
duction due to progenitor metallicity. In Figure 1131 we 
take the fiducial model and vary the abundance of stable 
^''Fe throughout the ^^Ni zone from zero to 25%, corre- 
sponding to a metallicity variation of roughly zero to 3 
times solar. The total mass in the ^^Ni zone is held fixed 
at 0.6 Mq - i.e., in this construction stable iron is pro- 
duced at the expense of ^^Ni. Thus the higher metallicity 
model has lower bolometric luminosity, lower tempera- 
tures, and an earlier onset of the 2 — » 1 recombination 
wave. The secondary maximum thus occurs about seven 
days earlier in the higher metallicity model. 

Because the size of the iron core remains fixed in the 
models of Figure lTSl one expects the absolute luminosity 
of the secondary maximum to remains nearly constant. 
This is indeed the case in the I and J band lightcurves. 
Interestingly, however, the H and K band secondary max- 
imum are slightly brighter in the low metallicity models. 
This is due to the significantly larger emissivity of cobalt 
compared to iron at these wavelengths (see Figure . 
Thus the relative strength of the secondary maxima in 
different NIR bands may turn out to be a useful measure 
of the ratio of cobalt to iron in the bulk ejecta, and hence 
the progenitor metallicity. 

Significantly, the effect of stable iron group elements 
on the NIR lightcurves is distinct from that of ^^Ni. The 
predicted correlation (i.e., earlier secondary maxima are 
as bright or brighter than later ones) conflicts with the 
primary observed trend. Variations in iron group ele- 
ments is thus one potential source of deviation from the 
standard NIR lightcurve behavior. This emphasizes the 
value of NIR observations in providing constraints of the 
progenitor metallicity and the early explosion processes. 

6.4. Calcium Distribution and the IR Triplet Feature 

On the whole, the distribution of intermediate mass 
elements in the ejecta has little direct effect on the NIR 
lightcurves. Emission from iron group lines dominates 
the continuum flux, and few individual IME line features 
are significant. The one exception is in the I-band, where 
the prominent Ca II IR triplet feature ((^/-weighted rest 
wavelength 8579 A) can potentially affect the photome- 
try significantly. The triplet lines are intrinsically very 
strong, and may remain optically thick far out in the low 
density layers of ejecta, even when the calcium abun- 
dance is very low (Xca « 0.01). 

Unfortunately, the net effect of the Ca II lines on the I- 
band lightcurve depends upon the details of the radiative 
transfer, specifically the line source functions. Two lim- 
iting cases can be considered: (1) If the calcium lines are 
assumed to be pure scattering, flux is nearly conserved 
over the line profile - i.e., the absorption component of 
the P-Cygni profile nearly cancels the emission compo- 
nent. In this case, the line feature should have only a 
minor effect on the broadband magnitude. (2) If, on the 
other hand, the calcium lines are assumed to be com- 
pletely absorbing (as in pure LTE calculations) the net 
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Fig. 14. — Effect of the Ca II IR triplet feature on the I-band 
lightcurve. For reference, the red dashed line shows the fiducial 
model lightcurve when all calcium is removed. If the calcium lines 
are assumed to be pure scattering, the IR triplet feature leads to a 
modest decrease in the I-band magnitude (solid black line). If the 
calcium lines are assumed to be pure absorbing, emission in the 
triplet feature increases the I-band magnitude significantly (blue 
dot-dashed line). 



emission in the Ca II IR triplet feature can be consid- 
erable, because radiation absorbed in blue part of the 
spectrum is allowed to re-emerge in the triplet lines. In 
the models of this paper we have always assumed pure 
scattering calcium lines, for the reasons described below. 

Figure ^1 quantifies the net effect of calcium on the 
I-band lightcurve by comparing the fiducial model to a 
similar model in which all calcium has been removed. 
Under the pure scattering assumption, the Ca II triplet 
feature causes modest (~ 0.1 mag) I-band differences. 
These arise for two reasons. First, because of blending 
with other lines, flux is not strictly cons erved o v er the 
triplet line profile. Second, because the iBesselll 1)19901) 
I-band filter profile cuts off near 8500 A, the absorption 
and emission component of the triplet feature are not 
equally sampled. Thus, under pure scattering, the cal- 
cium feature serves to slightly decrease the I-band mag- 
nitude. 

If, on the other hand, the calcium lines are assumed 
to be pure absorbing, net emission in the triplet features 
enhances the peak I-band magnitude by ~ 0.5 mag, and 
delays the time of peak by nearly seven days. The sec- 
ondary maximum is no longer distinct and prominent, 
and the model I-band lightcurve does not fit the observed 
lightcurve shape or colors. The result emphasizes that 
the I-band secondary maximum does not occur because 
of emission in the Ca II triplet feature, but rather in spite 
of it. 

From an atomic physics perspective, the low collisonal 
rates in the SN ejecta and the large oscillator strength of 
the Ca II triplet hues suggest that pure scattering is in- 
deed the more plausible representation of the line source 
functions. This is supported by our lightcurve fits to the 
SN 2001cl observations. To the extent that pure scat- 
tering does hold, the shape and peak brightness of the 
I-band lightcurves of SNe la will be generally insensitive 
(at the < 0.1 mag level) to the details of calcium produc- 



FlG. 15. — Dispersion in peak magnitude (measured at the first 
lightcurve maximum) as a function of wavelength band for the 
models of Figure [Tni with ^^Ni masses between 0.4 and 0.9 Mq. 



tion and distribution. This is an important condition for 
the robust application of I-band lightcurves to cosmology 
studies. One interesting exception is worth noting: if the 
calcium distribution deviates from spherical symmetry, 
flux scattered in the triple feature will be redistributed 
anisotropically along different viewing angles. 

The results of this section necessitate the inclusion of 
a non-LTE treatment of the calcium level populations 
within the fully time-dependent multi-dimensional radia- 
tive transfer calculation. This solution is in fact readily 
incorporated into the Monte Carlo approach, and will be 
included in the upcoming version of SEDQNA. The com- 
mon adoption of purely absorbing lines in previous LTE 
radiative transfer studies may be another reason why 
those calculations often found difficulty in fitting the I- 
band lightcurves of SNe la. 

7. NEAR-IR SN lA AS STANDARD CANDLES 

The cosmological interest in NIR observations of 
SNe la stems from the small observed dispersion (~ 
0.2 mag) in the peak JHK-band magnitudes. We ex- 
plore this property here using the models discussed in 
>i6.2l (and shown in Figure [TUIl for which the mass of ^^Ni 
is varied. We consider only those models with ^^Ni mass 
in the "normal" range of 0.4 to 0.9 Mq. 

Figure [T^ shows the spread in model peak magnitude 
(measured at the first maximum) as a function of broad- 
band filter. One notices a general decline in the dis- 
persion towards the red. In the B-band, the spread 
is roughly ±0.5 mag, similar to the o bserved (uncor- 
rected) dispersion l)Hamuv et al.lll99'6a|l . This decreases 
to ±0.4 mag in V and ±0.3 mag in R and I. In the NIR, 
the dispersion is more than a factor two in magnitude 
smaller than the B-band, being ^ 0.2 mag in J and K 
and ^0.1 mag in H. For the H and K bands, bear in 
mind the failure of the models to fix the observed first 
maximum, presumed an inadequacy of the atomic data. 
Nevertheless, this level of NIR dispersion is very similar 
to that inferred from observations (Mciklc 2000). 

In addition, the model peak NIR magnitudes show 
an impressive insensitivity to the other physical param- 
eters explored in this paper. Significant variations in 
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the progenitor metallicity, electron capture elements, and 
IME production all only led to small 0.1 mag) varia- 
tions in the absolute magnitude of the initial peak (Fig- 
ures ^1 ^1 and I14() . The mixing of ^^Ni also caused 
only slight magnitude variations unless extreme levels of 
mixing were considered (Figure Ej). These results further 
testify to the robustness of SNe la NIR lightcurves as 
standard candles. 

Our model lightcurves may also suggest further re- 
finements to cosmology studies in the NIR. The J-band 
lightcurves, for example, show almost zero intrinsic dis- 
persion at the local minimum occurring between the first 
and secondary maximum (about 15 days after B-band 
maximum; see Figure I1U|I . Even when the subluminous 
models are considered the magnitude variations at this 
epoch are small. It would be interesting to check if ob- 
served SNefe, exhibit a similar behavior. As it turns out, 
iMeikld l)2000() chose to study the J-band magnitudes of 
SNe la measur ed 14 days after B -maximum (as derived 
from fits to the'Eli as et alJ ljl985f) templates). He found 
that the dispersion in Mi4(J) was small, although not 
necessarily minimal. 

The low peak magnitude dispersion in the NIR fol- 
lows from the dependence of the model colors on the 
bolometric luminosity. The dimmer SNe also have lower 
temperatures, and radiate a greater percentage of the 
energy at redder wavelengths. This acts as a regulating 
mechanism which maintains a nearly constant peak mag- 
nitude in the NIR bands, regardless of the ^^Ni mass. On 
the other hand, the dispersion in ultraviolet and ?7-band 
magnitudes is reciprocally intensified, exceeding that of 
the bolometric lightcurves. 

8. DISCUSSION AND CONCLUSION 

We have modeled the far-red and NIR lightcurves 
of SNe la and given a detailed explanation of the 
characteristic secondary maximum. Our synthetic 
model lightcurves were calculated using parameter- 
ized 1-D ejecta configurations and the time-dependent 
multi-group radiative transfer code SEDONA. The model 
lightcurves displayed distinct and conspicuous secondary 
maxima and provided favorable fits to the NIR observa- 
tions of the normal Type la SN 2001el. By varying the 
mass of ^^Ni, the models also reproduced the observed 
trend that brighter SNe la have later and more promi- 
nent secondary maxima. 

We trace the origin of the secondary maximum directly 
to the ionization evolution of iron group elements in the 
ejecta. Specifically, the NIR emissivity of iron/cobalt 
gas peaks sharply at a temperature T21 « 7000 K, mark- 
ing the transition between the singly and doubly ionized 
states. The recombination of iron-rich gas from 2 —^ 1 
thus leads to enhanced redistribution of radiation from 
blue to NIR. Interestingly, as the supernova cools, the 
global ionization evolution takes the form a 2 ^ 1 "re- 
combination wave" gradually receding deeper into the 
ejecta. The onset and propagation of this wave through 
the iron rich layers marks the rise and fall the secondary 
maximum. Because the ejecta are transparent at longer 
wavelengths during these epochs, the NIR observations 
allow us to watch directly as the recombination wave 
scans progressively through deeper and deeper layers of 
ejecta. 

While the models considered in this paper captured 



the essential aspects of SN la NIR lightcurves, they also 
highlighted several outstanding issues for the radiative 
transfer calculations. First, the use of a complete and 
accurate atomic linelist proved critical in modeling the 
NIR bands. Further improvement of the currently avail- 
able atomic data is likely needed to accurately model 
the H and K-band lightcurves. Second, non-thermal ion- 
ization effects from radioactive gamma-rays become sig- 
nificant for igxp ?J 70 days when LTE predicts neutral- 
ity, and thus likely have a dramatic impact on the NIR 
lightcurves at these late epochs. Third, a proper treat- 
ment of Ca II IR triplet line source function is crucial 
in synthesizing accurate I-band model lightcurves, as the 
assumption of purely absorbing lines leads to unrealistic 
results. Future advances to the SEDDNA code will per- 
mit more detailed studies of these and other important 
effects. However, one does not expect the technical de- 
velopments to change the general NIR lightcurve trends 
and dependencies explained in this paper. 

In this paper, we studied the dependence of the NIR 
lightcurves on a number of important physical param- 
eters, which highlighted the many ways in which NIR 
lightcurves offer valuable diagnostics of the ejecta prop- 
erties and powerful constraints on explosion models. 

First, the double-peaked morphology of the NIR 
lightcurves can be taken as a direct consequence of the 
abundance stratification in SNe la, in particular the 
concentration of iron group elements in the central re- 
gions. This confirms previous inference s based upon 
post-maximum and nebular spectra (c.s.. Br anch et alJ 
1985; Hoflich et al. 2002; Kozma et al. 2005). Abun- 
dance stratification is generic to certain classes of explo- 
sion model s, for example the delayed detonation models 
llKhokhlo v II99ID and the detonation from failed deto- 
nation (^Plewa et alJ 1200*^ . In contrast, models char- 
acterized by large-scale mixing, s uch as published 3- 
dimensional defl agration models l|Gamezo et alJ 120031: 
iReinecke et al.l2002|') are likely inconsistent with the dou- 
ble peaked behavior in the NIR lightcurves. Further NIR 
observations, coupled with the transfer models, should be 
useful in constraining the exact degree of ^^Ni mixing in 
SNe la, and thus in testing current and future explosion 
paradigms. 

Second, the luminosity of the secondary maximum pro- 
vides a measure of the amount of iron group elements 
(both stable and radioactive) synthesized in the explo- 
sion. The observed correlation between B-band decline 
rate and the luminosity of the secondary maximum pro- 
vides strong evidence that slower declining SNe la have 
(on average) a larger production of iron group elements. 
This conclusi on coincides with several other inferences 
to the same iC ontardo et al.ll2000l iMazzah et "all 119981: 
iStritzinger etah 2006) . 

Third, the NIR secondary maximum is sensitive to 
the amount of stable iron group elements produced 
in the explosion, an d hence the progenitor metallicity. 
iTimmes et al.l 1)20031) has suggested that metallicity vari- 
ations may lead to 25% variations in ^^Ni. We find that 
the NIR signature of this variation is distinct from that 
of varying ^^Ni independently. Our predicted correla- 
tion arising from metallicity variations is that earlier sec- 
ondary maxima will be as bright or brighter than later 
ones. This conflicts with the primary observed trend, 
and suggests that the metallicity is likely a sub-dominant 
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cause of SN la luminosity variations. This signature fur- 
ther opens the possibihty of using an observational sam- 
ple of NIR lightcurves to constrain metallicity variations 
among different progenitor populations. 

Fourth, the timing of the secondary maximum is a di- 
rect probe of the temperature and ionization evolution 
in the ejecta, and hence a useful tool in interpreting ob- 
servations. As it turns out, the ionization evolution of 
crucial importance to the optical lightcurve decline rate 
as well. The 2 — > 1 recombination of iron group elements 
contributes significantly to blanketing in the B and V- 
bands due to the development of blends of strong Fe II 
and Co II lines in the optical part of the spectrum. In 
this context, NIR observations may also be especially 
useful in diagnosing peculiar objects. For example, the 
paradoxical SN 2002cx was ~ 2 mag subluminous, but 
had a peculiar spectru m characteristic of "hot" iron-rich 
ejecta l|Li et al.ll2003ft . NIR observations would compli- 
ment that spectroscopic analyses be providing indepen- 
dent measures of the temperature scale, ionization evo- 
lution, iron core size and mixing. 

Finally, we reiterate the cosmological utility of the NIR 
lightcurves of SNe la. In our lightcurve models, the dis- 
persion in NIR peak magnitudes is found to be small 
(< 0.2 mag) even when the physical properties of the 
ejecta are varied widely. At certain epochs (e.g., the 
J-band local minimum occurring near day 15 after B- 
maximum) the predicted dispersion is even smaller. Our 
models thus further testify to the robustness of SNe la 



NIR lightcurves as standard candles. 

Perhaps more importantly, the NIR observations pro- 
vide deep and readily interpreted diagnostics of the phys- 
ical conditions in the SN la ejecta, and therefore offer an 
opportunity for controlling potential systematic errors in 
the cosmology studies. We have shown that the physical 
parameters affecting the NIR secondary maximum are 
also those that may modulate the optical peak brightness 
and decline rate of SNe la. Thus the NIR lightcurves are 
one of the most promising places to search for empiri- 
cal secondary parameters related to deviations from the 
standard width-luminosity relation. 

In this context, the NIR observations strike a nice bal- 
ance in their level of information content. Given the 
double-peaked morphology, the NIR lightcurves carry 
considerably more information than do the single-peaked 
optical band observations; at the same time, they re- 
main simple enough to submit to a much more ready 
statistical analysis than do the complex spectral time se- 
ries. For all these reasons, a sizable sample of observed 
NIR lightcurves would be a remarkably powerful tool in 
constraining the possible evolution or progenitor drift of 
SNe la arising in differing stellar populations. This po- 
tential should be exploited further with expanded obser- 
vational programs and concerted theoretical efforts. 



The author thanks Kevin Krisciunas and Peter Meikle 
for helpful comments and suggestions. 
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